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We discuss the prospects for detection of high energy neutrinos from dark matter annihilation 
at the Galactic centre. Despite the large uncertainties associated with our poor knowledge of the 
distribution of dark matter in the innermost regions of the Galaxy, we determine an upper limit on 
the neutrino flux by requiring that the associated gamma-ray emission does not exceed the observed 
flux. We conclude that if dark matter is made of neutralinos, a neutrino flux from dark matter 
annihilations at the GC will not be observable by Antares. Conversely, the positive detection of 
such a flux would either require an alternative explanation, in terms of astrophysical processes, or 
the adoption of other dark matter candidates, disfavouring the case for neutralinos. 



PACS numbers: ... 

I. INTRODUCTION 

There is robust observational evidence for the domi- 
nance of non-baryonic dark matter over baryonic matter 
in the universe. Such evidence comes from many inde- 
pendent observations over different length scales. The 
most stringent constraint on the abundance of dark mat- 
ter comes from the analysis of CMB anisotropics. In par- 
ticular, the WMAP experiment restricts the abundance 
of matter to lie in the range iluh? = 0.135i";o°g Ij. The 
same type of analysis constrains the amount of baryonic 
matter to be in the range ^bh"^ = 0.0224 ± 0.0009, in 
good agreement with predictions from Big Bang nucle- 
osynthesis 0.018 < VLbh'^ < 0.023 (e.g. Ref. 0). 

It is commonly believed that such a non-baryonic com- 
ponent could consist of new, as yet undiscovered, parti- 
cles, usually referred to as WIMPs (Weakly Interacting 
Massive Particles). It is intriguing that some extensions 
of the standard model of particle physics predict the exis- 
tence of particles that would be excellent DM candidates. 
In particular great attention has been recently devoted 
to candidates arising in supersymmetric theories. The 
lightest supersymmetric particle (LSP), which in most 
supersymmetric scenarios is the so-called neutralino, is 
stable in theories with conservation of R-parity, and can 
have masses and cross sections of typical WIMPs. 

One possible way of probing the nature of dark matter 
particles is to look for their annihilation signal Q. For 
this purpose, the best regions to examine are those where 
the dark matter accumulates, the annihilation rate being 
proportional to the square of the particle number den- 
sity. A wide literature exists discussing the prospects of 
observing annihilation radiation from the Galactic cen- 
tre (e.g. Refs. 0, |Mflli,high energy neutrinos from 
the Sun (e.g. Refs. 01 B gamma-rays and syn- 
chrotron from dark matter clumps in the galactic halo 
(e.g. Refs. |lol ITll ll^ ll^) , gamma-rays from external 
galaxies (e.g. Refs. Il4l Il5l lia il7J 1. positrons and an- 
tiproton (e.g. Refs. [ll^nll20n2l|l and more. 

Large uncertainties are associated with predictions of 
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annihilation fluxes, due to our poor knowledge of the dis- 
tribution of dark matter, especially in the innermost re- 
gions of the Galaxy. Numerical simulations suggest that 
the dark matter density is well approximated by "cuspy" 
profiles, with a power-law behaviour p oc r~'^ . Estimates 
of 7 vary between having no cusp, 7 0, j22| |. to a cusp 
7=1 that is further steepened by adiabatic compres- 
sion of the baryons ^23] . One can trace these differences 
in large part to uncertainties in the stellar mass in the 
inner galaxy as inferred from microlensing experiments. 
The poor knowledge of 7 implies uncertainties of several 
orders of magnitude in the annihilation flux. The situa- 
tion is made even worse by the possible influence on the 
dark matter profile of the probable adiabatic formation of 
the supermassive black hole lying at the Galactic centre. 
Such uncertainties make indirect searches less effective 
for constraining the physical parameters (such as mass 
and cross sections) of dark matter particles. 

We suggest here a method for evading the astrophys- 
ical uncertainties in the neutrino flux, by requiring that 
the associated gamma-ray emission does not exceed the 
flux observed by the EGRET experiment in the direc- 
tion of the Galactic centre. In fact, if we normalize the 
gamma-ray flux to the EGRET data, the corresponding 
neutrino flux will be an upper limit on the actual neu- 
trino flux measurable on Earth. Choosing the EGRET 
normalization corresponds to fixing the product JavN^, 
where the quantity J, defined below, includes all of the as- 
trophysical information, av is the total annihilation cross 
section and is the number of photons produced per 
annihilation. 

This paper is organised as follows: we first discuss the 
gamma-ray source observed by the EGRET satellite in 
the direction of the Galactic centre ; in Sec. Ill we briefly 
review the results on the distribution of dark matter from 
observations and N-body simulations. In Sec. IV we 
present the particle physics details of our candidate, the 
neutralino, arising in supersymmetric theories, in Sec. V 
we review the prospects of indirect detection of such can- 
didates through gamma-ray and neutrino emission, for a 
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typical dark matter profile, and in Sec. VI we compare 
the prospects of indirect detection through annihilation 
radiation from the GC with other searches. We present 
in Sec. VII the upper limit on the neutrino flux, obtained 
by normalizing the annihilation flux to the EGRET data, 
and wc finally give our conclusions in Sec. VIII. 

II. THE EGRET SOURCE AT THE GALACTIC 
CENTRE 

The Galactic centre region has been observed by 
EGRET, the Energetic Gamma Ray Experiment Tele- 
scope, launched on the Compton Gamma Ray Observa- 
tory in 1991, and sensitive to an energy range 30MeV- 
30GeV. A strong excess of emission was observed in an 
error circle of 0.2 degree radius including the position 
I = 0°, b = 0°, the strongest emission maximum lying 
within 15 degrees from the GC 24]. 

The radiation exceeds, and also is harder than, the 
expected gamma ray emission due to the interaction of 
primary cosmic rays with the interstellar medium (see 
e.g. Strong et al. 1998 25]). At the energies we are inter- 
ested in, > 1 GcV, the main source of photons is the 
decay of tt*' mesons originating from processes such as 

P + X ^ TT° 



He + X ^ tt" 

where X is an interstellar atom. The interested reader 
will find a detailed estimate of the background radiation 
in Cesarini et al. 2003 26]. 

It is intriguing to conjecture that such excess emis- 
sion could originate from dark matter annihilation at the 
Galactic centre . However, such an interpretation is prob- 
lematic. In fact, as noticed by Hooper and Dingus [2^ . 
the EGRET source is not exactly coincident with the 
Galactic centre , which would make the interpretation of 
the signal as due to the annihilation in a spike around 
the Galactic centre at least problematic. 

Furthermore there is some evidence, although weak, 
that the source could be variable. Such a result could 
rule out completely the interpretation of the excess emis- 
sion as due to annihilation radiation from the Galactic 
centre . The variability of 3EG J1746-2851 has been re- 
cently discussed in Nolan et al. 2003 l2q. An additional 
flaw has been pointed out by P.Salati [5^, namely the 
fact that the HI column density was merely interpolated 
in the region of interest, where it was thought to be un- 
reliable due to strong self-absorption and high optical 
thickness. It is an open question how the conclusions 
would change if different assumptions are made about 
the HI column density. 

Here, we will regard the EGRET observation as an 
upper limit on the annihilation gamma-ray flux from the 
Galactic centre. 



HI. DARK MATTER DISTRIBUTION 

The usual parametrization for dark matter density pro- 
files is 



'^^ ^ (r/i?)-r[l + (r/i?)"](/'-'r)/" ' ' 

where r is the galacto-centric coordinate, i? is a charac- 
teristic length and a, (3 and 7 are free parameters. 

There is consensus, at present, about the shape of the 
profile in the outer parts of halos, but not in the inner- 
most regions, due to loss of numerical resolution in N- 
body simulations and to the poor resolution in observa- 
tion of rotation curves of outer galaxes. Navarro, Frenk 

6 White [23, found with N-body simulations that the 
profile could be well approximated at small radii with a 
power-law p(r) = with 7 ~ 1. Othe r gr oups reached 
different conclusions (see e.g. Refs. [30ll3]]]) . 

The most recent N-body simulations [33 . l3^ Is^ sug- 
gest that profiles do not approach power laws with a well- 
defined index at very small radii. Profiles continue to be- 
come shallower, i.e. the (negative) logarithmic slope be- 
comes higher, when moving towards the centre. Some au- 
thors however contend that convergence is reached with 

7 « 0.2 - 0.3 at 0.3% of the virial radius [13 . 

An additional complication of the dark matter pro- 
file at the center of our Galaxy, is the well-established 
presence of a a 3.6 x 10^ solar mass black hole (see 
e.g.Ref. 113]), that would accrete dark matter, producing 
a so-called 'spike' and leading to an enhancement 
of the annihilation fiux by several orders of magnitude 
(see Ref. |i6|j, and references therein, for a discussion of 
indirect detection of dark matter in presence of spikes, 
and of dynamical effects that could potentially destroy 
them). 

The observational situation is even less clear. The 
analysis of rotation curves of galaxies has led some au- 
thors to claim inconsistency of the observed 'flat' pro- 
files with the cuspy profiles predicted by N-body simu- 
lations. Other groups 0, ^3 claim instead that cuspy 
profiles are compatible with observations. Hayashi et al 
|32| compared the observational data directly with their 
numerical simulations (rather than fits of their simula- 
tions) and found no significant discrepancy in most cases. 
They attributed the remaining discrepancies to the dif- 
ference between circular velocities and gas rotation speed 
in realistic triaxial halos. 

It is clear that the predictions of annihilation fiuxes are 
strongly affected by the uncertainties in dark matter dis- 
tribution. In particular, since the annihilation rate is pro- 
portional to the square of the particle density, different 
profiles can lead to uncertainties of many orders of mag- 
nitude. To get around these, we will use the gamma-ray 
flux observed by EGRET in the direction of the Galac- 
tic centre , to get rid of astrophysical uncertainties and 
produce a robust upper limit on the neutrino flux from 
dark matter annihilation at the Galactic centre . 
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IV. SUPERSYMMETRIC DM: NEUTRALINOS 

Neutralinos are by far the best studied dark mat- 
ter candidates. They arise in supersymmetric theories 
with conservation of R-parity, in which the hghtest su- 
persymmetric particle (LSP) cannot decay in standard 
model particles, and is thus stable. In most cases the 
LSP is the neutralino, i.e. a linear combination of the 
supersymmetric partners of the gauge and higgs bosons 

X(= x?) = znS + z,2W3 + z^^H^ + z.^Hl (2) 

The matrix z diagonalizes the neutralino mass matrix, 
which is expressed as 

Ml —mzcpsw mzspsw 

M2 mzCjSCw —mzSfjCw 

—mzcpsw mzCfjCw — /i 

mzspsw -mzspcw -fJ- 

(3) 

in the basis {B, W3, H^, H^). 

Similarly defining Vii(2) as the wino (higgsino) frac- 
tion of the lightest chargino, the neutralino annihilation 
channels and cross-sections most relevant for indirect de- 
tection are 
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Annihilation in these channels thus increases with the 
wino or higgsino fraction of the neutralino. The spec- 
tra of the indirect detection signals studied here keep an 
imprint of the dominant channel. 

For muon via neutrino production: the W~^W~ and 
Zh channels produce more energetic neutrinos, i.e a 
harder neutrino spectrum than bb. Both the neutrino- 
nucleon cross section ((Tiz-at) and the muon range {Rp.) 
being proportional to neutrino energy, harder spectra 
give higher muon detection rates for the threshold con- 
sidered here (5 GeV): (j)^ oc <j)^(T^^M{E^)Rf^[E„). 

For gamma production: the bb and also tt channels 
dominate the 7 spectra around 2 GeV but at higher en- 
ergies, the harder WW and Zh channels come in. Ex- 
periments with different thresholds can thus see different 
processes. 

The influence of the dominant annihilation channel is 
displayed on figure below. 

We have performed a scan of SUSY models at the GUT 
scale, computing renormalisation group equations and 
radiative electroweak symmetry breaking with Suspect 
[36j, the neutralino relic density with Mi cr omegas 
and detection rates with Darksusy 38| 59] . The SUSY 
models explored fall in 3 classes (see j8| for definitions) 

CMSSM: with universal scalar mg and gaugino mi/ 2 
mass parameters in the ranges: 

50GeV < mo < 4000GeV, 50GeV < mi/2 < 2000GeV, 
^0 = 0, tan/3 = 5,20,35 



Non universal gaugino mass M2\gut' same values as 
above, except for A'^Igut — 0.6mi/2 (instead of ITO1/2), 
leading to AI2 ~ Mi in the neutralino mass matrix (eq. 
|3Jl; the resulting non-zero wino contents (212) allows for 
non-negligible relic densities. 

Non universal gaugino mass M^Igut^ same values 
as in the universal case (tan /3 — 20, 35 only) with 
MsIgut = 0.6toi/2 (instead of lmi/2); to decrease the fi 
parameter in the neutralino mass matrix (eq. ISJ to favour 
the higgsino fraction (^13(4)) and decrease scalar masses, 
in particular the pseudo scalar A mass. We do not re- 
lax Higgs sector universality, whose interesting effects on 
dark matter are similar to those of a lower M^Igut- 

Finally, we apply the following conservative cuts on 
our models: 

• Higgs mass: rrih > 113.5 GeV [3^ . 

• Chargino mass: > 103.5 GeV [iol ]. 

• Rehc density: 0.03 < fi^h^ < 0.3, but we also show 
the WMAP Q range fig^o^^^h^ ^ 0.1126^°;°}^}, 

• 6 ^ S7 Constraint 41]: 

2.33 X 10-4 < BR{b sj) < 4.15 x lO"", 

• The muon anomalous magnetic moment |56| : 
8.1 10-1" < (5™"y = Sl^'P - (5™ < 44.1 10-1" [2a]. 

Given the recent evolution of this last range, the on- 
going debate about the use r-decay data [53 and the 
drastic effect of this 2a cut, which both excludes the 
SM and many interesting dark matter models, the range 
< S^'^'^'y < 8.1 10"i" will not be discarded, but displayed 
in pale on all plots. 



V. GAMMA-RAY AND NEUTRINO FLUX 
FROM THE GC 

Indirect detection of Dark Matter is based on observa- 
tion of annihilation products like gamma-rays, neutrinos 
or synchrotron emission of secondary electron-positron 
pairs. The spectrum of secondary particles of species 
i from annihilation of DM particles whose distribution 
follows a profile p{r) where r is the Galacto-centric co- 
ordinate, is given by 



dN, 



1 



dE AttM^ Jiine of sight 



dsp'^ (?-(s,-0)) 



(5) 

where the coordinate s runs along the line of sight, in a 
direction making an angle ■0 respect to the direction of 
the GC. av and dNi/dE are respectively the annihilation 
cross section and the spectrum of secondary particles per 
annihilation, while M is the mass of the annihilating DM 
particle. 

To isolate the factor depending on astrophysics, i.e. 
the integral of along the line of sight, we introduce. 
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following the quantity J{ip) 



.5 kpc V 0.3 GeV/cm^ / Jlme of sight 



dsp'^ ir{s,ip)) 
(6) 



and its average over a spherical region of solid angle AQ, 
centered on ip = 0, J {Ail). 

With these definitions the flux from a solid angle AH. 

is 



<^>^{A^l,E) ~ 5.6 X 10 



-12 



m 

dE 



M 

X Aficm^^s 



(7) 



Apart from astrophysics, large uncertainties on the 
quantities in eq.Hare associated with the details of parti- 
cle physics. The dependence of the annihilation cross sec- 
tion on the mass is different for each DM candidate, 
and even in the framework of a specific supersymmetric 
scenario, cross sections for a given mass could span over 
several orders of magnitude. 



and consider the EGRET source as an upper limit on 
the WIMP annihilation flux. In this sense, we see from 
Fig. n that if neutralinos are the dark matter particle, 
then there is room for profiles even more "cuspy" than 
NEW. 

Always assuming a NEW profile, we show in Eig.|21the 
neutrino-induced muon flux from dark matter annihila- 
tion at the GC. We show for comparison the expected 
sensitivity of the Antares telescope (e.g. 48], currently 
under construction in the Mediterranean sea. The tele- 
scope sensitivity depend on the incoming neutrino spec- 
trum, we thus show two sensitivity curves (for a 3-years 
period of observation), one relative to a hard flux (rele- 
vant for the W'^W~' and Zh channels), the other relative 
to a soft flux (relevant for the bb channel). As can be seen, 
the predictions fall several orders of magnitude below the 
Antares sensitivity. Of course, at this stage, this does not 
necessarily imply that Antares will not observe any neu- 
trinos from the Galactic centre , as we have seen in the 
previous section that it is possible that the actual dark 
matter profile is steeper than NEW, adopted for Eigs. 
and la 
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FIG. 1: Gamma-ray flux from neutralino annihilation at the 
GC, assuming a NFW profile. For comparison we show the 
EGRET and GLAST sensitivities. Shades paler than in the 
legend denote a low S^^'^^ value. 

We show in Fig.^the gamma-ray flux from neutralino 
annihilation at the Galactic centre assuming a NEW pro- 
file, along with EGRET and GLAST sensitivities. We see 
that all the supersymmetric models predict fluxes below 
the EGRET sensitivity in this case, but many of them 
could produce fluxes observable by GLAST. 

Nevertheless, as already mentioned, EGRET did ob- 
serve a source at the Galactic centre , although it is un- 
clear whether this emission is actually to be attributed to 
WIMP annihilations. We adopt a conservative approach 
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FIG. 2: Neutrino-induced muon flux from neutralino annihi- 
lation at the GC, assuming a NFW proflle. For comparison 
we show the expected Antares sensitivity. Shades paler than 
in the legend denote a low 5^"^^ value. 



VI. COMPARISON WITH OTHER SEARCHES 

In this section we compare, for completeness, the 
prospects of detection of the SUSY models discussed 
above with other detection techniques, which are actually 
insensitive to the profile of dark matter in the innermost 
regions of the Galaxy. 
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In Fig. Owe show the flux of neutrinos from neutraHno 
annihilation in the solar core. The projected sensitivities 
of both Antares and IceCube appear to be able to 
probe the supersymmetric models with the non-negligible 
higgsino fraction necessary for an efficient neutralino cap- 
ture rate in the Sun. 
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900 



FIG. 3: Neutrino-induced muon flux from neutralino anni- 
hilation in the solar core. Shades paler than in the legend 
denote a low 5^"*^^ value. 



We also show in Fig. 0| the potential of direct detec- 
tion techniques to probe the neutralino nature through 
the search for neutralino-nucleon interactions in large de- 
tectors, such as Edelweiss 50] and CDMS 51j. There 
are a couple of orders of magnitude between the present 
-day experiment sensitivities and the most optimistic pre- 
dictions for neutralinos. But this gap could be bridged 
by next-generation experiments such as Edelweiss II (e.g. 
lil) and Zeplin 



VII. UPPER LIMIT FOR THE NEUTRINO 
FLUX 

In order to maximize the neutrino flux from dark mat- 
ter annihilation at the Galactic centre , we normalize the 
flux of gamma-rays, associated with such a neutrino flux, 
to the EGRET data. This corresponds to fixing, for each 
model, the product JavN^, with — ^^'^re 
Ri is the branching ratio of all the channels i contributing 
Ni gamma-rays above a given threshold energy. 

Having fixed the particle physics contents of our dark 
matter candidate, the ratio between the number of pho- 
tons and the number of neutrinos emitted per annihila- 
tion is known. We can thus estimate the neutrino flux 
from the Galactic centre associated with a gamma-ray 
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FIG. 4: Prospects of neutralino direct detection. Shades paler 
than in the legend denote a low 5^"^'' value. 



emission reproducing the EGRET data. Finally we can 
convert the flux of neutrinos into a flux of muons, pro- 
duced by neutrinos interactions with the rock around de- 
tectors on Earth, in order to compare with experimental 
sensitivities. 

The rescaled flux of muons (f>^^°™{> Eth) will thus be 
given by 



REGRET (j;^ 



(8) 



where the label NFW reminds that NFW profiles have 
been used to compute profile- independent flux ratios, and 
-B* is the energy at which we decide to normalize the flux 
to the gamma-ray data (in our case E^, = 2GeV). 

The results are shown in Figs. [S] and |3 The muon 
flux normalised to the EGRET data represent an upper 
limit, as the observed gamma-ray emission could be due 
to processes other than dark matter annihilation. The 
comparison with the Antares sensitivity shows that only 
the highest mass neutralinos can possibly be detected in 
the Galactic centre. Insisting on the WMAP relic den- 
sity in Fig. Eland using the hard neutrino spectrum sen- 
sitivity appropriate to the relevant Zh channel, we need 
at least 700 GeV neutralinos, whose contribution to the 
muon anomalous moment is similar to the (excluded?) 
Standard Model. 

If neutrinos are nevertheless observed above the given 
fluxes, then their interpretation as due to neutralino an- 
nihilation is problematic and would actually require ei- 
ther the adoption of other dark matter candidates an- 
nihilating dominantly into neutrino pairs or a different 
explanation, e.g. in terms of astrophysical sources. 
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FIG. 5: Neutrino-induced muon flux from the Galactic cen- 
tre normalized to EGRET with relic density values. Shades 
paler than in the legend denote a low SJ^^^ value. 



FIG. 6: Neutrino-induced muon flux from the Galactic cen- 
tre normalized to EGRET, for models with WMAP-preferred 
relic density sorted by leading (= BR > 0.5) annihilation 
channel. Shades paler than in the legend denote a low S^^^ 
value. 



Concerning other dark matter candidates, a case- by- 
case analysis is needed. For Kaluza-Klein candidates 
(e.g. Ref. and references therein) , in particular, 
there are several channels contributing to the neutrino 
flux (see Ref. H^). Neutrinos coming from the decay of 
charged pions originating in quark fragmentations have 
a relatively soft spectrum, and cannot be detected with 
Antares, even normalizing the gamma-ray flux to the 
EGRET data. A similar conclusion applies for neutri- 
nos from prompt semi-leptonic decay of secondary heavy 
quarks, despite the fact that the spectrum in this case is 
harder. One last channel could be potentially interest- 
ing, the direct production of neutrinos, which is nearly 
forbidden in the case of neutralinos. This channel is par- 
ticularly interesting since in this case the spectrum of 
neutrinos is a line, at energy equal to the mass of the 
Kaluza-Klein particle. Rescaling the fluxes obtained in 
Ref. we estimate this flux to be comparable with the 
Antares sensitivity to line spectra. A detailed analysis of 
this case will be presented elsewhere. 

Finally, to show how our upper bound on the neutrino 
flux from the Galactic centre would evolve with new data 
on gamma ray fluxes, we show in Fig. the flux above 
60 GeV coming from the same neutralino annihilations in 
the Galactic centre, applying the same normalization to 
EGRET that we used for neutrinos in Fig. [S] As in that 
figure, the points trace an upper bound on the gamma 
flux above 60 GeV, given the EGRET measurement. If 
Hess sees a signal (which is not excluded according to 
Fig.Q), e.g. two orders of magnitudes below this gamma 
upper bound, the upper bound on the neutrino flux Fig.[Sl 
can accordingly be reduced by two orders of magnitudes. 



VIII. CONCLUSIONS 



The flux of neutrinos from dark matter annihilation at 
the Galactic centre depends on the assumed dark mat- 
ter profile and on the details of annihilation of the specific 
candidate adopted. It is nevertheless possible to obtain 
an upper limit for the neutrino flux, by requiring that 
the associated gamma-ray emission do not exceed the 
flux observed by EGRET in the direction of the Galactic 
centre . 

We have estimated such upper limits in the case of 
neutralinos and concluded that any associated neutrino 
flux lies below the experimental sensitivity of Antares, 
unless the neutralino mass is above ^ 700 GeV . In this 
case, corresponding to models with a low 6^^^^ value, and 
even assuming that the gamma-ray emission observed by 
EGRET is entirely due to neutralino annihilation, the 
upper limit on the neutrino flux is barely above the min- 
imum signal observable by Antares in 3 years. 

This means that Antares will not be able to see neutri- 
nos from neutralino annihilation at the Galactic centre . 
Conversely, the positive detection of such a flux would 
either require a different explanation in terms, e.g., of 
other astrophysical sources, or the adoption of dark mat- 
ter candidates other than neutralinos. 
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